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ABSTRACT 

The spin probability distribution of Dark Matter haloes has often been modelled as 
being very near to a lognormal. Most of the theoretical attempts to explain its origin 
and evolution invoke some hypotheses concerning the influence of tidal interactions 
or merging on haloes. Here we apply a very general statistical theorem introduced 
by Cramer (1936) to study the origin of the deviations from the reference lognormal 
shape: we find that these deviations originate from correlations between two quantities 
entering the definition of spin, namely the ratio J/M^^^ (which depends only on mass) 
and the modulus E of the total (gravitational + kinetic) energy. 
To reach this conclusion, we have made usage of the results deduced from two high 
spatial- and mass resolution simulations. Our simulations cover a relatively small vol- 
ume and produce a sample of more than 16,000 gravitationally bound haloes, each 
traced by at least 300 particles. We verify that our results are stable to different sys- 
tematics, by comparing our results with those derived by the GIF2 and by a more 
recent simulation performed by Maccio et al. 

We find that the spin probability distribution function shows systematic deviations 
from a lognormal, at all redshifts z < I. These deviations depend on mass and redshift: 
at small masses they change little with redshift, and also the best lognormal fits are 
more stable. The J — M relationship is well described by a power law of exponent a 
very near to the linear theory prediction {a = 5/3), but systematically lower than this 
at z< 0.3. We argue that the fact that deviations from a lognormal PDF are present 
only for high-spin haloes could point to a role of large-scale tidal fields in the evolution 
of the spin PDF. 

Key words: methods: N-body simulations - methods: numerical 



1 INTRODUCTION 

Within the hierarchical galaxy formation model. Dark Mat- 
ter (hereafter DM) haloes are thought to play the role of 
gravitational building blocks, within which baryonic dif- 
fuse matter collapses and becomes detectable. On galactic 
scales, the formation of stars and their evolution provides 
an important p robe of the evolution of the visible conten t 
of the Universe (| White fc Reesll 19781 : IWhite fc Freii^llQQlf ). 
although the subtleties of the stellar formation processes 
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within galaxies, as of today not yet completely understood, 
hinders an exploitation of these objects as a clean probe 
of the evolution of DM haloes. On the other extreme of 
the mass scale, the most massive clusters of galaxies are 
regarded as one of the most reliable cosmological probes 
(|Bahcall et al.lll997^ : in particular, their abundance and evo- 
lution with redshift is a v ery sensitive test of the under- 
lying cosmological model (I White et all 1199 3': 'Hcnr vlll997l: 
Bahcall et all Il997l: [Oukbir fc BlanchardI Il997: Eke et al. 



199S : iDonahue fc Void Il999l: iHenrvl I2OO0I : lArhipova et al 
2003 : lHenrvll2004l : iGladders et al.ll2007l ) . 

The Mass Function (hereafter MF) is an indirect test of 
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the total virialised mass of DM haloes: exact predictions 
of the latter can be done using t he nonlinear spherical col- 
lapse model (jGunn fc Gotdl 19721 ). an essential ingredient of 
the Press-Schechter model for the MF. However, DM haloes 
also possess angular momentum, and taking into account 
its effect on the gravitational collapse of rotating shells has 
been shown to have a detectable consequence on the MF 
(|Del Popoldl2006l . l2009h . Thus, an investigation of the dis- 
tribution of angular momentum, and its connection with 
mass is particularly useful, as high-resolution simulations of 
increasing resolution will produce MFs with a very small 
statistical uncertainty. 

An exact determination of the shape of the DM halo spin 
PDF can also have important consequences for the abun- 
dance of low surface brightness galaxies, if t he latter form 
preferen tially within hig h-spin DM haloes jjimenez et al.l 
[l998 : Maccio et al.|[2007l ). Finally, it is of great importance 
also in models of the formation and evolution of central 
Black Holes. In collapse models where most of the baryons' 
specific angular momentu m is a fixed fraction of that of their 
host dark matter haloes (|Fall fc Efstathioulll980l : IMo et^ 
1 19981 ). the angular momentum and extent of the gaseous 
central accretion disc ar e strongly dependent on their's spin. 
IVolonteri fc Ree^ (|2005h find that the central density of the 
disc varies as po — thus the initial rate of accretion 

of the central Black H ole turns out to be a very sensitive 
function of the spin A (iPeebleslll97ll ). 

For all these reasons, investigations of the origin of the an- 
gular momentum growth and of the spin PDF of DM haloes, 
and of their evolution, can have an impact on many different 
open issues in large-scale str ucture format ion and evolution. 
Theoretical investigations (|Peeblesl [l969l ) predict that to 
zeroth-order the angular momentum should have a power- 
law dependency on the total virialized mass, with exponent 
5/3. They also make a prediction concerning the spin, a di- 
mensionless quantity defined as: 



They alternatively suggest an empirical fit: 



(1) 



where J =\ J \, E =\ Eki,z + Egrav \ are respectively the 
moduli of the angular momentum, and of the total (ki- 
netic plus potential) energy. The PDF of A has been pre- 
dicted to have an approximately lognormal distribution. 
This initial prediction was subsequently found to be basi- 
cally val id also when higher orde r effects were taken into 
account (jCatelan fc Theunslll996al lbh . More recently, the ra- 
dial p rofile of A has been d erived from modified Jeans' equa- 
tions (|Schmidt et al.|[2009l ). and found to be in good agree- 
ment with results from simulations. 

The properties of the spin distribution of DM haloes have re- 
cently been considered, taking advantage of the availability 
of hig h spatial- and mass-resolution simulations. iBett et al] 
l|2007h have analysed a large sample of haloes drawn from 
the Millennium simulation, and found (among other things) 
that the global spin distribution is poorly described by a 
lognormal distribution: 
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which has the shape of a lognormal except at very high and 
very low values of the spin. They also suggest that the actual 
shape of the distribution depends on the adopted numerical 
definition of halo. 

Theoretical and numerical studies are aimed at understand- 
ing the origin of the (almost) lognormal spin distribution, 
using a limited set of statistical and dynamical assump- 
tions or by perform ing controlled numerical experiments. 
IChiueh et al.l (|2002l l have shown that the distribution of spe- 
cific angular momentum (which they define as j' = J/M^^^) 
can be described by a rather complex PDF, which can be 
approximated by a lognormal in the central part, but de- 
viates significantly from it at low- and high values of j' . 
They also make predictions for the dependence of the peak 
of the spin distribution for different values of mass, suggest- 
ing t hat it varies linearly wi th the average crs of the distribu- 
tion. [Knebe^^fc^^oro^ (|2008l ) find a correlation of spin A with 
mass, albeit a weak one. Note that they use a set of simula- 
tions using boxes of varying sizes, while in the present work 
we adopt a single simulatio n, thus minimising the spurious 
tidal effects. Furthe rmore, (jBett et all l2007l : ICao fc: Whit3 
I-Mlll1 :| Li et al.ll2008l ). studying the spin distributions at low 
redshifts, have found that more massive haloes show larger 
values of A. Extensio ns of these results to high er redshift 
have been provided bv lDavis fc NataraianI (|2009l '). who have 
studied the evolution of the spin distributions at 2 > 6, and 
showed that more massive haloes (M ~ W^Mq) tend to 
have a median A higher than that of M ~ IO^Mq haloes. 
Also, high-A haloes tend to cluster more (by a factor 3- 
5) than low-spin haloes, a trend which strengthens with 
time. However, their simulation is restricted to a small box 
[Lb = 2.46/i^^Mpc), thus making their result more prone to 
uncertainties from cosmic variance. 

The steady improvement of the available hardware and soft- 
ware resources makes today possible simulations where the 
limits imposed by the finite spatial and mass resolution lim- 
its are challenged. The simulations we have performed in this 
work were aimed at providing a statistically significant sam- 
ple of reasonably well-resolved DM haloes. We have obtained 
a catalogue of more than 77600 DM haloes, each resolved by 
more than 20 particles, and we have used only those haloes 
with more than 300 particles, thus resulting in a catalogue 
containing more than 16400 haloes. We have chosen a box 
size of L = 70 Mpc, smaller th an the one used in pre- 
vious papers (e.g. IShaw et al.ll2006l ). and a large number of 
particles, to maximize the mass resolution. 
In this work we present the results of a numerical experi- 
ment, where we study the evolution of the angular momen- 
tum/spin distributions within a relatively small cosmologi- 
cal volume. We find that the spin distribution shows some 
clear deviations from a lognormal, and does not seem to re- 
lax towards the latter, above statistical uncertainty. These 
deviations do not seem to depend significantly on the redshift 
or the mass of the halo, at least for the underlying LCDM 
model studied in this paper. We also perform a possible 
check on the role of possible systematics in our analysis by 
compa ring with the result s of two different simulations: the 
GIF2 (|Giocoli et al.|[2"008h . and a simulation recently per- 
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fomed by iMaccio et alj (|2008l ). Both these works adopted 
a Spherical Overdensity (SO) halo finder, slightly diflerent 
from our chosen halo fi nder (the AmiRa Halo Finde r, AHF), 
recently introduced bv lKnollmann fc Knebd (|2009l V We ad- 
dress the problem of the origin of these deviations by ex- 
ploiting an_exact result from statistics {Cramer's theorem, 
ICramer Il936l ) to demonstrate that the deviations from a log- 
normal distribution are induced from correlations between 
the total energy and mass of haloes, as should be expected 
for not completely virialised haloes. 

The plan of the paper is as follows: In section [2] we present 
the details of the simulations we have performed, briefly de- 
scribing the code adopted, the initial conditions and the halo 
finder. We then describe the main results concerning J — M 
relationship (section [3TT}. In section|3]we concentrate on the 
spin probability distribution, and we show that the devia- 
tions from a lognormal shape are dependent on halo mass. 
We trace the origin of these deviations to the presence of de- 
tectable correlations among the quantities entering its defi- 
nition, by applying the Kendall test. We discuss our findings 
in the Conclusion (section [S} . 

In the following, we will denote the natural logarithm us- 
ing the symbol "In" , and the decimal logarithm as "log" (or 
"Log"). 



2 NUMERICAL SIMULATIONS 

As underlying cosmological model we ha ve chosen the 
fiduc ial 5-year WMAP LCDM cosmology (|Dunklev et all 
|2009| ). with Hubble constant: Ho = 71.9, and the 
cosmological parameters: (f2(, , Q,a , f2m , n^ , as) = 
[0.044, 0.742, 0.258, 0.963, 0.7986]. 

Our main target was that of producing a statistically sig- 
nificant number of DM haloes, thus we tried to adopt a 
softening length sufficiently small to ensure that the small- 
est haloes we want to resolve have a size large enough at the 
initial redshift of the simulation. In order to do this, we de- 
fine at first an initial instability radius riniuai , for which we 
refer t o the r2oo radius presented in eq. (1) of lMo fc White! 
(|2002D : rc-nn = 2.98x10'^ {Mg/Q^y^'-^ / (1+z) kpc (where Mg 
is defined as Mq). For the initial redshift of our sim- 

ulations, their formula gives: rinuiai = '"20o(-z = 50, Mm) ~ 
10.36 kpc. Moreover, as we explain later, we consider only 
haloes having a lower mass threshold Mm > 1.44 x IO^Mq. 
We then choose a softening length / = 12.5 kpc, thus allow- 
ing the linear density field to become gravitationally unsta- 
ble on this spatial scale at the beginning of the runs. Note 
that our use ri„itiai is circumscribed only to the initial in- 
stability criterion. 



2.1 Initial conditions and code 

We have chosen a box with size Lt = 70 Mpc, not large 
enough to minimise cosmic variance effects. Within this box, 
we have performed two runs, differing only for the number of 
particles used: run 32M was performed using 320'^ particles, 
and run 500M with 800^ Thus, the particles masses for runs 
32M and 500M are, respectively, 7.49 x 10* and 4.79 x 10^ 
h.-^ Mq. 

We have generated the initial conditions using a parallel 



version of the IC package by ISirkd (|2005l ). that we devel- 
oped ourselves. We started all the simulations from a redshift 
z — 50. The latter is chosen in such a way to ens ure that the 
linea r modes are linear within the chosen box (|Lukic et al.l 
|2007| . sec. 4.1). Recent work (|Knebe et al] l2009h suggests 
however that, for the mass range of interest to this work, 
the starting redshift has little infiuence on the average final 
properties of the haloes. 

Our main simulation tool is FLY, a parallel MPI N- 
body cosmological simulation code implementing a par- 
allel version of the Barnes-Hut octal tree algorithm 



(|Antonuccio-Delogu et "alll2003l : iBecciani et"al]|2007f ). FLY 

adopts an efficient parallelization scheme: domain decom- 
position is applied to distribute particles, and workload de- 
composition is furth er applied to ensure load balancing (see 
IBecciani et ahllioOOl . for a detailed description of the parallel 
algorithm implemented in FLY). 

2.2 Extraction of haloes 

The availability of large samples of DM haloes, obtained 
from state-of-the art N-body simulations, means we can ob- 
tain statistics such as the mass function with very small 
poissonian errors. However, in order to extract a significant 
amount of cosmological information from these statistics, 
one should have a control on the systematic effects. The 
definition of DM halo is one possible source of systematics. 
Presently, halo definition algorithms fall into two broad cat- 
egories: those based on the friends of friends (FOF) algo- 
rithm, and those based on some threshold overdensity cri- 
terion, both of which are usually supplemented by some 
recursive scheme to eliminate outliers (i.e. gravitationally 
unbound particles). Different definitions mostly affect small 
haloes and the outer, low -density regions of more massive 
haloes (|Lukic et al.ll2007l ). For this reason, we have chosen 
to restrict our attention to haloes described by a relatively 
large (300) minimum number of particles. 
We have decided to adopt for our analysis a group finder 
which has only one threshold paramet er, i.e. the re- 
cently introduced AMIGA Halo Finder (|GiU et al.ll2004l : 
iKnollmann fc Knebd [2OO9I ). In AHF haloes are defined us- 
ing the commonly used virialization criterion, i.e. first com- 
puting isodensity contours, and then including only those 
particles whose average density is larger than the critical 
oversdensity: 



Pc{rv)/Pb = Ayir{z) 



(4) 



In the above equation pc is the average density within the 
virial radius of each halo, pi, the background density, and 
the virial overdensity Avir{z) depends on the cosmological 
model. This is not the only possible physicall y meaningful 
definition of the mass of a numerical halo fsee IWhit3l200ll . 
for a discussion): the most relevant differences in these def- 
initions arise in the identification of the outer boundary of 
a halo, which can in principle have a significant impact also 
on the derived statistical properties. As we stated above, we 
minimize discreteness effects by including only haloes having 
a number of particles A'j, ^ 300, thus restricting the mass 
range to M > 2.247 x lO" and 1.437 x 10^" Mq, for the 32M 
and 500M runs, respectively. 

In AHF the gravitational energy of each particle is com- 
puted using the same tree data structure used in the octal 
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Figure 1. Evolution of the J — M relationship. The continuous 
lines are the best fits for two redshifts: z = 13.7 {upper) and 
z = 0.3 {lower), and the open symbols represent averages for 
the binned intervals. Error bars are poissonian. The mass of the 
haloes is evaluated as M„ir- 



tree codes to compute the gravitational force. Thus, the po- 
tential energy of haloes is computed without making any 
hypothesis concerning the shape or density distribution of 
the haloes themselves. 



3 ANGULAR MOMENTUM DISTRIBUTION 

Due to the rather limited size of the box we have employed, 
cosmic variance could have an influence on the statistical 
properties of the halo sample. In this section we will then 
first investigate the angular momentum-mass relation. 



3.1 Angular momentum— mass relation 

One of the most important predictions of tidal torque the- 
ory is that the angular momentum-mass relationship of re- 
laxed haloes sh ould be describ ed by a power law: J oc M", 
with oi = 5/3 (|Peebles|[l969l ). We test this prediction on 
the 500M run, shown in Figure [T] We notice that the distri- 
bution has an apparent excess of haloes having an angular 
momentum larger than the best-fit value. This is more evi- 
dent a.t z = 13.7, where the fraction of unrelaxed, high-spin 
haloes is larger, but we notice that also ai z — 0.3 the trend 
seems to persist. 

The best-fit values for the power law coefficients are shown 
in Figure [2] and we notice immediately that they are statis- 
tically compatible with the theoretical value, although the 
value at z = 0.3 is smaller (at the la level). There is then a 
deviation from linear theory predictions at late times, when 
structure growth was non-linear. One possible explanation 
could be that linear theory does not take into account merg- 
ers, an d their role in determin ing the statistical evolution of 
haloes (|Vitvitska et al.l bOOd ). In the same Figure we also 
show the effect of changing the definition of halo mass on 
the fitted properties of the J — M relationship. While in the 
lower part, and in the rest of the paper, we use the virial 



Figure 2. Best fitted exponents of the J-M relationship. The up- 
per plot show s results obta i ned u sing the static halo definition 
introduced bv lCuesta et al ] l|2008l '). the lower one is for the stan- 
dard AHF halo mass definition. The horizontal continu ous line 
show the theoretical a = 5/3 v alue from linear theory ijPeeblesI 
1 19691 : ICatelan fc Theunslfigoeal lbl) . Error bars show the standard 
deviations in the values of the fitted slopes. Although the upper 
values are slightly larger, the difi'crcnccs arc within the uncertain- 
ties. 



mass as deduced by AHF, in the upper part of this figure we 
show the effect of using a different definition of halo mass, 
recently introduced by Cuesta et al. (2008). As one can ap- 
preciate, the effect is very little. 



4 DISSECTING THE SPIN PDF 

The evolution of the angular momentum distribution of 
dark matter hal oes ha s often been modelled using the tida l 
torque theory jHoyld Il95ll: ISciamal 1 19551: IPeeblesI Il969l : 



lDoroshkevichlll970l : IWhitelll984l : Tschafedl2009h . which pre- 
dicts that a typical DM halo gains most of its angular 
mome ntum during the linear regime (|Catelan fc TheunsI 
Il996al ): it has further been shown that contributions from 
later collapse epoch s and non-linear corrections are small 
(jChiueh et al.ll2002l ). 

The s pin parameter A (eq. [T]) was introduced by iPeebled 
(|l97lf ) as a convenient dimensionless quantity to character- 
ize the tidal growth of the angular momentum, also through 
the study of its probability distribution function (PDF), 
P(A). All the quantities entering the d efinition of A dep end 
ultimately on the adopted halo finder (|Bett et al.ll2007l ). 
In order to understand the origin of this deviation, we will 
first consider the evolution of the spin PDF. In Figs [3]- [5] 
we show the evolution of the PDF for different mass inter- 
vals. We also plot the best-fit lognormal distributions, and 
in Table [T] we present the parameters of thes e fits. Although 
the fitting functions of lBett et all (|2007f ) and lHiotelid (|2008l ) 
provide a better fit than a lognormal, the latter is a phys- 
ically motivated PDF distribution for the origin of angu- 
lar momentum, predicted both by tidal torque and merger 
models. Thus, we will try to understand the possible phys- 
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Figure 3. Spin PDF at z=l, for different mass intervals. Error 
bars are Poisso nian, computed according to the prescriptions in 
iHeinrichI | |2003D . The curves are best-fitting lognormal functions 
(see Table [T] for the values of A and ui^gx)- 



ical origin of the deviations from the lognormal distribu- 
tion. It is apparent from these figures that there is deficit of 
high- A haloes (log A/ A > 1.3 -r- 2o- i ogA)i at all redshifts, as 
was observed pr eviously ('Gardnei| [200ll : IShaw et al.l l2006l : 
iBett et al.ll2007l : [Hiotolis, ,200a ') . This feature does not sig- 
nificantly depend on the mass, and in the next section we 
will argue that it originates from statistical correlations be- 
tween quantities entering the definition of spin, and from 
deviations from lognormal behaviour of the energy E. 
From Table [T] we also notice that the location of the max- 
imum A depends on the mass, being a slowly decreasing 
function of mass. This fact has been previously noticed 
by various authors (ICole fc Laceylll996l : iMaccio et al.ll2007l : 
iBett et al.ll2007l : iKnebe fc PowedbOOSl ). but they suggested 
that the A— M relationship fiattens for z < 1, while our sim- 
ulations show instead that also at more recent epochs the 
relationship holds true. Our fi ndin gs and those of Kricbc and 
Power are different to those of lLem sori & Kauffmann (1999,), 
but more consistent with those o f|M accio et all (j2007f l. who 
find a very mild decreasing trend. Finally, we notice an evo- 
lution also in the shape of the best-fitting lognormal func- 
tion, in the sense that the dispersion aiogx seems to increase 
with decreasing redshift (third column of Table [ij . We do 
not have sufficient statistics for z > 0.1 to determine if in 
the highest mass bin the high-A deficit and the aiogx width 
follow the same trend. 

We summarize all these features in Figure[6l which shows an 
interesting fact: in the high mass bin 

the location of the maximum A stays almost constant after 
z = l. 

Although these deviations from the lognormal fit are be- 
yond statistical uncertainty, one could still suspect that they 



Figure 4. Same as Fig.O z=0.5. The PDF for the highest mass 
bin is not plotted, due to the poor statistics. 



Table 1. Parameters of the lognormal fits of the PDFs shown in 
Fig-EJE] Columns are as follows: (1) Mass interval (in Mq/i"^ ), 
(2) A (3) X (4) redshift 
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could be a product of some systematic effects, for instance 
introduced by the adopted halo finder. We have then checked 
the spin P DF distributions for two different simulations: the 
GIF2 run (IGiocoli et al ] 120081) ■ and a more recent simula- 
tion performed by Maccio et al.l (|2008l ) . We only use the 
halo catalogues provided by these simulations, instead of 
directly applying the AHF to their raw data. In the GIF2 
haloes were identified using a spherical overdensity criterion, 
without pruning the halo of gravitationally unbound parti- 
cles. Also Maccio et al. use a spherical overdensity criterion, 
but they select haloes based on a number threshold criterion 
Np ^ 250. 

While Giocoli et al. use the same definition of spin as we do 
(e g. [H, Maccio et al. use instead the defintion introduced 
bv lBullock et al.l (120011 ). and we convert it to ours using the 
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Figure 5. Same as Fig.[3l for redshift z=0.1. 
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Figure 6. Evolution of the lognormal fits to the spin PDF, for 
different mass bins. The different curves show the fits for different 
redshifts: continous: z = 0.1, dotted: z = 0.5, dashed: 2 = 1. 

prescirption described after eq. 

In Fig. [7] we show the PDFs for these simulations. The most 
striking feature of this figure is surely the great discrepancy 
in the distribution width (rms o"iogA) of run 500M wrt the 
other runs. Differences in the quantity A are less pronounced, 
though still present. The values we find are however consis- 



Figure 7. Comparison of spin PDF at redshift z=0.1 for dif- 
ferent simulations. The continuous curves are best-fitting log- 
normal fun ctions for three diff erent runs: run 50 0M ( continu- 
ous, fe/acfcl. lMaccio et al.l | |2008|) (dotted, red), and lGiocoli et al] 
l|2008l ^ (dashed, blue). 




-3 -2 -1 



Figure 8. Same as Fig. [T] but here the dotted red line uses only 
the relaxed haloes from the simulation of Maccio et al. are plotted. 
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tent with those previously fou nd by other autho rs, as re- 
cently summarised in Fig. 7 of IShaw et alj (|2006l ) (see also 
their section 4.1). The reasons of such a result are still under 
investigation by us. Perhaps, the cause of these differences 
is to be found in the different ratio between the amount 
of accretion- or merger-do minated spin histories (but see 
iD'Onghia fc NavarrolTioOTl ) . Another interesting feature is 
that the PDFs of both Giocoli et al. and Maccio et al. show 
an excess of low spin haloes, and a deficit in the high-A tail 
with respect to a lognormal, while the PDF for the 500M run 
only shows the latter featu re. A similar exc ess of low-spin 
haloes was observed e.g. bv lBett et all (|200it ). and seems to 
be a general feature of many other simulations. Note how- 
ever that the mass range of our simulations is more biased 
towards galaxy-sized hal oes, while low-spin halo es tend to be 
those with higher mass (|Knebe fc Poweij|2008h . The excess 
of low-spin haloes seems to be increasing with the average 
mass range probed by the simulation, which is larger in the 
GIF2 simulation, and smaller in our run 500M. On aver- 
age, low-spin, high-mass haloes could be less relaxed at any 
given epoch: this could explain why our haloes seem to bet- 
ter follow a lognormal distribution. 



iMaccio et al.l (|2008l ) introduce a sophisticated criterion to 
isolate relaxed haloes (see section 2.2. of their paper): how- 
ever, we see from Fig. [8]that these relaxed haloes have a spin 
PDF showing the s ame qualitative fe atures as the global one. 
Also, we note that lBett et al.l (|2007h found that picking just 
relaxed haloes made a big difference on the form of the spin 
distribution, but that the biggest difference was due to the 
halo-finder rather than imposing a cut on the instantaneous 
"virial ratio". Thus, the difference between relaxed and un- 
relaxed haloes could be very slight, in the light of these dif- 
ferences. 

We then conclude that the spin PDF at all redshifts since 
z ~ 0.5 shows a detectable deficit of haloes having log A/A > 
1.3a"iogA, with respect to a lognormal distribution. This 
deficit is quite detectable and statistically robust, and ap- 
parently is a feature not depending on the particular halo 
finder or dynamical state of the halo. 

In the next subsection we will use a very general statistical 
theorem to understand the origin of this feature in the spin 
PDF. 



4.1 Statistical origin of deviations 

The spin A can be regarded as a stochastic variable depend- 
ing on three other variables: J, E and M. We can decompose 
the spin in two terms, plus a constant: 



ln(A) = In 



J 



ln(G) 



(5) 



The first factor in the right-hand side can be regarded as de- 
pending on the mass M, in the light of the scaling J cx M'^^^, 
while the total energy E also depends on the dynamical 
state of the halo. The latter is also determined by envi- 
ronment, through e.g. tidal fields and streaming motions. 
We expect these two terms to be distributed similarly for 
fully virialised haloes, although the environmental depen- 
dence of E could introduce an additional dependence in ad- 
dition to that on the mass M. Here however we only con- 
sider the distribution of the global spin, without analysing 
environmental differences, but the reader should be made 




5 6 7 

ln(J/M»(5/2)) 




ln(j;M«(5/2)) 

Figure 9. Correlation diagrams for the quantities in the right- 
hand side of eqn. [S] at rcdshift z= 1. Here and in Figures 10-12 
axes units are arbitrary. The four levels are density eontours as 
fraetions of the maximum, at the 80%, 60%, 40% and 20%, levels, 
from the inner- to the outermost. The three plots are for different 
mass intervals, corresponding to those of the PDF distribution 
of Fig. [5] upper left: 10 ^ log(M) < 10.75, upper right: 10.75 < 
log(M) s: 11.50, lower left: 11.50 < log(M) < 12.25. 




Figure 10. Same as Fig.|9l for z=0.5. 



aware of the fact that environmental differences are not only 
modulate d by the mass dep endence, as recently shown by 

[Faltcnba cher fc White! (|201(]| '). 

Following Cramer's theorem (|Crameilll93^ ). two necessary 
and sufficient conditions for the variable ln(A) to be nor- 
mally distributed are that the two quantities in the right- 
hand side of eq.[5]are statistically independent and normally 
distributed. Thus, we expect that deviations from a lognor- 
mal fit arise if at least one the two terms in eq. [S] is not nor- 
mally distributed, and/or if there are correlations between 
J/M^^^ and E. In our case, both situations are taking place. 
In order to check the independence of the two terms, we show 
in Figure[9t lll [ contour plots of the quantities J/M^^^ and E. 
It is evident that the two quantities show correlations in all 
mass bins, although these become progressively uncorrelated 
with increasing mass. The lack of correlation in the larger 
mass bin (bottom right plot in Fig. Illf) is probably partly 
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ln(j;M«(5/2)) 



ln(J/M"(5/2)) 



Figure 11. Same as Fig. [9] for z=0.1. In addition to the tliree 
mass intervals, in the , lower right section we plot the contour 
plots for the highest mass bin; 12.25 < log(M) ^ 13. 



Table 2. Kendall's correlation indexes tj^ between J/M^/^ and 
E. A value of the correlation index lying in the interval [—0.5, 0.5] 
is taken as evidence of poor correlation. 



logM ([Mq]) 


tk 


z 


10.00 - 


10.75 


-0.391 


1 


10.75 - 


11.50 


-0.400 


1 


11.50 - 


12.25 


-0.329 


1 


10.00 - 


10.75 


-0.382 


0.5 


10.75 - 


11.50 


-0.381 


0.5 


11.50 - 


12.25 


-0.427 


0.5 


10.00 - 


10.75 


-0.368 


0.1 


10.75 - 


11.50 


-0.388 


0.1 


11.50 - 


12.25 


-0.391 


0.1 


12.25 - 


13.00 


-0.297 


0.1 



due to the lack of very large haloes in our sample, a conse- 
quence of cosmic variance. The o utputs of Kendall's correla- 
tion test (|Kendall fc Stuartll 19791 ') are given in Table They 
show that the two quantities become progressively more sta- 
tistically independent with increasing mass, and that the 
correlations are actually never very significant, in statistical 
terms. This is also evident from the plots, which show that 
the two outermost contours (40 and 20% of the maximum), 
are partially elongated. Although small, this correlation can 
partially explain the presence of the deviations of the spin 
from a lognormal distribution. 

Finally, we have checked whether one or both of the 
two quanties entering our decomposition of the spin A, i.e. 
J/M^/'^ or E, are lognormally distributed. In Figure [T^ we 
show the PDF of log(i5), which shows a shoulder at rela- 
tively large values, at least up to M < 3.6 x 10^^ Mq. For 
larger masses, the distribution deviates even more, a sign 
of poor relaxation of these haloes. Note that the dynami- 
cal relaxation of haloes generally takes longer, affecting the 
gravitational binding energy and thus the total energy, so 
that this could explain the lack of a lognormal distribution 
of the latter. 
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Figure 12. PDF of the total energy E, for z= 1. (continuous), 
z= 0.5 (dotted), z= 0.1 (dashed). We do not show the PDF for 
the highest mass bin, which is noise dominated. 



5 CONCLUSIONS 

The work we have presented in this paper aimed at un- 
derstanding the origin of the deviations of the spin PDF 
from a lognormal distribution. Such deviations have been re- 
cently noticed in numerical simulations, and even predicted 
by statistical mode ls for the angular momentum growth 
(|Chiueh et al.ll2002l) . Instead of attempting some ad-hoc fit- 
ting function for the PDF, we have attempted to understand 
the statistical origin of these deviations. To that end, we 
have perfomed two high-resolution simulations (maximum 
resolution; mp = 4.37 x 10^ Mq), in a box of Li, = 70h 
Mpc, which provided us with a large number (~ 16, 000) 
of well resolved haloes, i.e. haloes having more than 300 
particles. Thus, we believe that our results are statistically 
robust. 

We can summarise our main findings in the following points: 

• The J-M relation is well fitted by a power law, with 
exponent a = 5/3, except at very recent epochs; 

• The spin PDF is systematically low at log A/ A > 

1.3 (TlogA- 

We have concentrated our attention on intepreting the lat- 
ter point, using basic statistical arguments (Cramer's theo- 
rem). An implicit assumption in our approach is that the 
lognormal distribution should be in some way preferred 
over other possible functional forms. Notice that not only 
the "neares t neighbours" tidal torque theory predicts such 
an output (jCatelan &c T hcuns"l9 96al Q), but als o merging- 
based models like those of Vit vitska et al.l (|2002l ) a nd, more 
recently, the extended Press-Schechter models by iHiotelid 
(|2008D . Some of these models were based on the analysis of 
N-body simulation outputs, with a mass resolution down 
to rup ~ 10^ Mq, but using a set of nested boxes (e.g 
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iMaccio et allbOOTl 'l. In order to probe halo formation us- 
ing a comparable mass resolution, while avoiding the tech- 
nical problems arising from the correct reproduction of tidal 
fields, that arise when using nested boxes, we have per- 
formed a larger simulation, and analysed the evolution of 
the mass f u nction . 

iBett et all (|2007l ) have provided arguments to support the 
idea that a PDF which is not anti-biased against small 
(A ~ 0) values of the spin (a typical feature of the log- 
normal distribution) should provide a better description of 
the actual spin PDF. However, our results do not seem to 
support their conclusion: particularly at low values of A the 
lognormal seems to provide a very good fit of our simula- 
tions. One possible reason of this discrepancy could be due 
to the lies in the fact that our simulations, compared to 
the Millennium run analysed by Bett et al., do not have 
haloes more massive than ~ IO^^Mq, while in the Millen- 
nium run this upper limit extends to M « 10^^ M0. If high 
mass haloes tend to have a smaller spin (|Knebe fc Power! 
I2OO8I ). then we conclude that our simulations probe better 
the mass range typical of low- to intermediate mass galactic 
haloes, avoiding instead massively galactic and cluster-sized 
haloes. The latter tend to be not completely relaxed, thus 
their spin could deviate from the lognormal form expected 
from Tidal Torque Theory. 

We have chosen a relatively small cosmological volume 
(Li, = 70/i~^Mpc), in order to test whether the angular 
momentum-mass relation and spin distributions show sta- 
tistically detectable deviations from the predicted analytical 
forms they have for dynamically relaxed haloes. The large 
statistics made possible by the high resolution allows us to 
detect some significant deviations: however, any numerical 
experiment is subject to some systematics. 
We cannot rule out that some residual systematic uncertain- 
ties could arise from the adopted halo finder: the actual spin 
of a numerical halo should also depend on the adopted halo 
finder, i.e. on the numerical definition. Even in the "n ear- 
est neighbours" tidal torque theory by iPeeblej (|l969l ) the 
distribution of angular momentum within the (spherically 
averaged) collapsed halo reflects the evolution of the torque 
during the collapse: the torques on shells placed at larger 
distances decrease on average with radius. Thus, one could 
expect that the distribution of A would depend on the ra- 
dius at which it is estimated, and this is generally different 
for halo flnders based on critical isodensity contours or on 
friends-of-friends algorithms. 

Concerning the environmental dependence of the spin dis- 
tribution, it is now e vident that the original claim by 
iLemson fc KauffmannI j 19991 ). who did not find any evi- 
dence for this dependence, was due to a lack of sufficent 
numerical re solution in their simulations. A lready in a pre- 
vious work I Antonuccio-Delogu et ahl |2002| ) we found that 
the spin's PDF has a significantly smaller dispersion and 
mean in Voids, a result consistent with the more recent 
findings that a higher merger frequency tends to produce 
haloes with higher spi n (lGardneill200ll : lMaller et al.ll2002l : 
iHetznecker fc Burkertlbooel '). 

Our simulations have sufficient statistics to test the predic- 
tions of the linear model concerning the J-M relationship. 
We notice that the index a of the power law J oc M" at red- 
shift zero is smaller than the standa rd value 5/3. A s imilar 
fact had been previously noticed bv lSugerman et al] (|2000l . 



see their Table 6), using simulations on a larger volume, thus 
it is not probably a peculiarity of our use of a small volume. 
We measure average values of the spin at 2 ~ 0.1 consistent 
with those found recently by other authors jMaccio et al.l 
[2003): A ~ 0.03. Maccio et al. find little dependence of A on 
halo mass, even when they restrict their analysis to relaxed 
haloes. Our results are c onsistent with their findings, and 
with theoretical analyses jChiueh et al.l [20021 '). We do how- 
ever see a dependence of the spin dispersion on mass: more 
massive haloes tend to have a lower dispersion, although it 
is a very slight trend which decreases towards low redshifts, 
consistently with the idea that it is induced by the presence 
of not fully relaxed haloes. 

The deviations of the spin's PDF from a lognormal arise 
because of statistical correlations between mass and angu- 
lar momentum and because of significant deviations of (at 
least) the total energy E from a lognormal distribution, as 
implied by Cramer's theorem. These correlations tend to 
disappear when haloes collapse, as one should expect for 
ideal completely relaxed haloesij. Note that this conclusion 
is true both for merger- and collapse-induced relaxation: in 
either case, the result is that the halo retains no "memory" 
of its previous state, due to dynamical friction acting on 
dynamical timescales. 
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